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ABSTRACT 

Recent high-quality observations of low surface brightness (LSB) galaxies have shown that 
their dark matter (DM) halos prefer flat central density profiles. On the other hand the standard 
cold dark matter model simulations predict a more cuspy behavior. One mechanism to reconcile 
the simulations with the observed data is the feedback from star formation, this might be suc- 
cessful in isolated dwarf galaxies but its success in LSB galaxies remains unclear. Additionally, 
including too much feedback in the simulations is a double-edged sword, in order to obtain a cored 
DM distribution from an initially cuspy one, the feedback recipes usually require to remove a 
large quantity of baryons from the center of galaxies, unfortunately they also produce twice more 
satellite galaxies of a given luminosity than what is observed. Therefore, one DM profile that 
produces cores naturally and that does not require large amounts of feedback would be preferable. 
We find both requirements to be satisfied in the scalar field dark matter model. Here, we consider 
that the dark matter is an auto-interacting real scalar field in a thermal bath at temperature T 
with an initial Z 2 symmetric potential, as the universe expands the temperature drops so that 
the Z2 symmetry is spontaneously broken and the field rolls down to a new minimum. We give 
an exact analytic solution to the Newtonian limit of this system and show both, that it satisfies 
the two desired requirements and that the rotation curve profile is not longer universal. 



Subject headings: galaxies:formation-galaxies:halos 
6946)-galaxies:fundamental parameters 

Introduction 



The longstanding core/cusp discussion, whether 
the central dark matter (DM) profiles in dwarfs 
and low surface brightness (LSB) galaxies are 
more core-like and rounder than the standard 
cold dark matter ( CDM) model predi c ts, re - 
mains an open iss ue (jvan Evmeren et al. (2009); 
see de Blok ( 2010f ) for a recent review). So far, 
the core profiles most frequently used in the lit- 
erature an d that best fi t the observations are 
empi rical ( Burkert 1995 ; iKuzio de Narav et al. 
2010h . Though they are useful to characterize 



-galaxies:individual (NGC 1003, NGC 1560, NGC 



the cores, especially since more and more recent 
high-quality observations of LSB galaxies suggest 
that the core-like behavior (p ~ r -°- 2 ) j s p re _ 
ferred in the central reg i ons of dwarf and LSB 
galaxies (lOh et al.l lioilt iRobles fc Matosl 12012 ; 



de Blok et al.l 120011 ). This is a problem to the 
CDM model which prefers p ~ r _1 at small r 



( Navarro et al. 2010() . Thou gh the latest s i mula- 
tions can reach p ~ r -°- 8 dNavarro et al. 201dl; 



properties of galaxies, it is necessary to find a 
theoretical framework which naturally produces 
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Merrit et allbOOri iGraham et al.ll2006l ). which is 
not in total agreement with observations. 

The current trend to solve the core/cups dis- 
crepancy in the CDM model is to include the dy- 
namics o f the baryonic component into DM sim- 
ulations iGovernato et al.ll201ol l2012t iMacio et all 
201ll;IStinson et al.ll201ll;lRomano-Daz et al.l200cl) . 



By including feedback from star formation in sim- 
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ulatio ns of field dwarf galaxies, iGovernato et al 
(|2012n managed to change an initially cuspy DM 
halo into a core-like halo. On the other hand, 



Sawala et al.l (|2012l ) find in their simulations twice 
as many satellites of a given luminosity around a 
Milky Way size host halo. They found this by 
considering, in addition to feedback, the effect 
that surrounding galaxies have on its host halo by 
means of tidal interation. This means, satellite 
galaxies are too luminous. They also find that the 
DM halos of satellites are more strongly affected 
than their stellar component by tidal interactions, 
and that, in the case that cores are a universal 
feature of dwarf galaxies (as suggested by the ob- 
servations), it is difficult to reconcile the observed 
high total mass-to-light ratios in dwarf spheroidal 
galaxies (dSphs) with those of their simulated ha- 
los under the CDM paradigm. This discrepancy 
might be solved with a different way of galaxy 
formation. 

An issue that remains unanswered is whether 
the same feedback recipes used in dwarf galax- 
ies _wo£k__as__wBn_Jn_^jSBs^^ As p ointed out 
in iKuzio de Narav fc Spekkens ( 2011 ) this seem 
unlikely, it is necessary to show that there 
is an accord between the content of gas in 
LSBs (usually low) and the l arge amount of it 



needed in CDM simulations ferook et al. 2012: 



Stinson et all 1 2 1 lL IGuo et alJ 
20 lie IKuzio de Narav fc Kaufman 



. More et al 
201ll) . We 

know LSBs are a larg e portion of the total galax 



2010t iMoi 
tannl l201lt 



ies that are observed ( e.g.. McGaugh et al. 1995), 



therefore, as far as there is not yet an agreement 
with LSBs and CDM simulations it worths explor- 
ing alternative models that provide a solution to 
the above disagreements. 

There are several models in the literature that 



(Avila- Reese et al. 


2001; Ccmbranos et al. 


2005 


Strigari et al. 2007 


; Spergel & Steinhardd 


2000 



ify CDM, like the warm dark matter (WDM) 
and self-interacting dark matter (SIDM) models, 
haven't been able to solve these discrepancies yet 
dNavarro et al.ll2010t Kuzio de Narav fc Spekkens 



201ll:IZavala et al.l2009HDave et al.l200lHYoshida et al 



2000). There are models that modify gravit y like 
f(R) the ories foe Felice fc Tsui ikawal lSoToh and 



MOND (iMilgroml l2010b ISandersI 120091) . but they 
are currently more at the effective-theory level 



rather than at a fundamental one. Nevertheless, 
there are some rotation curves fits of LSBs galax- 
ies in MOND models whose fits are almost perfect 
to the observed data (| Sanders! 20091 ). 



One model that has received much atten- 
tion is the scalar field dark matter (SFDM) 
model. It is our aim to show that in this model 
there is an scenario of galaxy formation (de- 
scribed in section 2) different from the stan- 
dard model used in CDM simulations and that 
naturally produces core density profiles, repro- 
duces rotation curves of large and small galax- 
ies on equal footing as MOND and empirical 
dark matter models do, but that doesn't require 
large amounts of feedback. In previous works it 
has been verified that the SFDM model repro- 
duces cosmological observations as well as CDM 



(Ror 


rfguez-Montova et al.l 12010: Suarez & Matos! 


2011 


Robles & Matos 


20121: M 


aeana et al.l 2012at 


Bohmcr & Harko 2007 


; Harko 


2011). 



This article is organized as follows, in section 
2 we describe the SFDM model to be analysed in 
this paper, in section 3 we give our results and 
section 4 is devoted to discussion and conclusions. 

2. SFDM model 

2.1. Previous work and unsolved issues 



The idea was first considered by Sin (199 4J) and 
i ndepe ndently introduced by iGuzman fc Matosl 
(|200(Dh . In the SFDM model the main hy- 
pothesis is that the dark matter is an auto- 
interacting real scalar field that condensates form- 
i ng Bose-Einstein C ondensate (BEC) "drops" 
([Magana et al.ll2012al) . Then we interpret these 



BEC drops as the halos of galaxies, such that its 
wave properties and the Heisenberg uncertainty 
principle stop the DM phase-space density from 
growing indefinitely, and thus, it avoids cuspy 
halos and reduces the number of small satellites 
(|Hu. Barkana fc Gruzinovll2000l) . 

In the SFDM model, the scenario of galaxy for- 
mation proposes that galactic halos form by con- 
densation of a scalar field (SF) with an ultra-light 
mass of the order of m ~ 10 _22 eV ( units where 
the speed of light c=l). From this mass it fol- 
lows that the critical temperature of condensation 
Tcrit ~ Ti~ 5 / 3 ^TeV is very high, thus, they form 
BEC drops very early in the universe. It has been 
proposed that these drops are the halos of galax- 



2 



ies (jMatos fe Urenal 120011) . i.e., that halos are gi- 
gantic clumps of SF. The DM halos can be de- 
scribed in the non-relativistic regime, where DM 
halos can be seen as a Newtonian gas. When the 
SF has self-interaction, we need to add a quartic 
term to the SF potential and in the Newtonian 
limit the equation o f state of the SF is that of a 
polyt rope of index 1 (jSuarez fc Matosll201lt Harko 
20 111 ) . On the other hand, big str uctures form like 
in the CDM model, by hiera rchy (jMatos fc Ureha 
20011 ISuarez fc Matoslboill ). thus, all predictions 
of the CDM model at big scales are reproduced by 
SFDM. 

Different issues of DM halos and the cosmo- 
logical evoluti on of the SFDM model have bee n 
investigated in Colpi et al.l (Il986 ): 
Matos fc Urenal <|200lh : IChavanis 



Gleiserl (119881) : 
(|201lh . Al- 



though the model is notably succesfull, there are 
at least two reasons to complement the model 
in addition to finding an explanation to the two 
discrepancies discussed in the introduction. 

One of them is found when we consider the fully 
condensed system at temperatur e T=0, the fits to 
rotat i on curves (RCs) o f LSB s ( Robles fc Matod 



2012; IBohmer fc Harkol 120071) show deviations 



from the observed data at large radii because DM 
halos are modeled only by considering the ground 
state (complete condensation), whose DM density 
and v elocity profile are given by (jBohmer fc Harkol 
20071) 



p°(r) 
V 2 (r) 



sin(7rr/-R) 

Po 7~— > 
■kt/R 

kxGpl (sax{Kr) 
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Kr 



(la) 
cos(Kr) ) (lb) 



where K = ir/R, p§ = p°(Q) is the central density, 
and the halo radius determined by p{R) = is 



R = TT 



m 

Gm 3 



(2) 



with h is Planck's constant divided by 2n, m the 
mass of the DM particle, G the gravitaional con- 
stant, and b the scattering length. The density 
profile depends on two fitting parameters, the cen- 
tral density and a length scale R. However, R 
depends only on fixed parameters, the interaction 
parameter A and the mass of the SF particle, which 
implies that it should not vary from galaxy to 
galaxy. However, when fitting rotation curves of 



galaxies, as R is a fitting p arameter it gets differ- 
ent values for each gala xy (jRobles fc Matosll2012 ; 
IBohmer fc Harkoll2007h . this is something contra- 
dictory and represents a problem for the simplest 
SFDM model. 

A second problem lies in the fast decrease of 
the velocity profile (equation (|lbj) ) after its max- 
imum value, such decrease causes a disagreement 
between the fits and the observed data in large 
galaxies because the data usually remain "flat " 



until the outermost regions (|Swaters et alj|2000h 



In addition to this, and for large galaxies, when 
we aim for the best fit to the velocity maximum, 
we obtain a worse fit in the outer regions, i.e. the 
better the fit to the velocity maximum in the RC 
the worse the fit becomes in the outer regions and 
viceversa. 

One a pproach to solve t he pr oblems was con- 
sider by Matos fc Urenal ( 2007 ). and later by 
lHarko fc Madarassvl ( 2011 ). It consisted in in- 
cluding fin ite temperature of dark ma tter in the 
DM halos. lHarko fc Madarassvl (|201ll) show that 
by including a small correction to the pure con- 
densed state, albeit in a different way from the one 
in this paper, it is possible to solve the first and 
partially solve the second discrepancy. The inclu- 
sion of temperature T, does mainly two things: 1) 
it makes the halo radius temperature-dependent, 
thus it's not fixed for all galaxies anymore and 2) it 
lifts the RC fit in the outer region and keeps it flat 
until the last value. However, effect 2) is neglible 
when the halo temperature is T<0.5TbeCj where 
Tbec is the critical temperature of Bosc Einstein 
condensation, if Tbec ^TeV, we would expect 
present halos to be well approximated by equation 
(llbp and hence, we will be unable to simultane- 
ously obtained a good fit to the RC maximum 
and the RC outer regions. Though this last prob- 
lem is not readily visible in galaxies with small 
radius (outer radius of <10 kpc), it is notable in 
large galaxies, therefore problem 2) is just par- 
tially solved. 

Thus, an alternative approach to solve these 
two SFDM discrepancies considers non-condensed 



=0, 



SF configurations at T 
tions in exited states (|Balakrishna et al 



SF configura- 
" 19981: 



Urena fc Bemall2010l: iBernal et alj|20ld) . These 



configurations fit RCs up to the last data point 
and can even reproduce the wiggles se en at large 
radii in high-resolution observations (jSinl [1994; 
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Colpi et al.l [1986). The problem found in this 
type of solution was that the configurations re- 
quired for the fits (4-5 exited states) were unsta- 
ble a nd decayed to the ground state in a short 
time ( Siddhartha fc Urenall2003h . Therefore, we 
will expect to see today that most DM halos are 
in their ground state, which means that the dis- 
agreement at large radii would remain. 

2.2. Finite Temperature SF scenario 

Motivated to solve all these issues we consider 
the following scenario which includes temperature 
of the DM and the exited states of the SF. 

The idea is that the dark matter is a scalar field, 
we consider a real SF with a repulsive interac- 
tion embedded in a thermal bath of temperature 
T, we also consider the finite temperature corrrec- 
tions up to one-loop in the perturbations. This is 
described by the p otential (jKolb fc Turnerlll987 ; 
DalfoTOeEaD[i999) 



1 rh 2 c 2 



X 



-$ 4 



A . 



-kiT 2 <$> 



TT 2 kiT 4 



V ^ 2 h 2 ' ' 4 ' ' 8'"^ " 90ft 2 c 2 

(3) 

for the case when fceT>> fnc 2 . Here fee is Boltz- 
mann's constant, A = X/(h 2 c 2 ) is the parameter 
describing the interaction, (i 2 := m 2 c 2 /h 2 is a pa- 
rameter, and T is the temperature of the thermal 
bath. The first term in V($>) relates to the mass 
term, the second to the repulsive self-interaction, 
the third to the interaction of the field with the 
thermal bath, and the last to thermal bath only. 

At some high enough temperature in the early 
universe, when the SF reaches the minimum of the 
potential ($«0) it behaves as dust, reproducing 
all the features of the CDM model at cosmological 
scales. In this period the SF interacts with the 
rest of the matter, eventually, when the temper- 
ature is sufficiently small the SF decouples from 
the interaction with the rest of the matter and fol- 
lows its own thermodynamic history, cooling down 
with the universe expansion. After this moment 
and for sufficiently low T, the term proportional 
to T 4 is not longer important and can be dropped 
out. Morover, as the temperature continues de- 
creasing, the initial minimum of the potential in 
$=0 eventually becomes a local maximum, after 
this happens the initial Zi symmetry of the po- 
tential V($) is broken. The latter happens at a 



critical temperature Tc given by 

2mc 2 



Vx 



(4) 



The critical temperature T^ determines the mo- 
ment in which the DM fluctuations can start grow- 
ing, they do it from the moment when T<Tc until 
they reach a stable equilibrium point, for example 
in $ 2 nin =fc 2 i (T^ - T 2 )/4 (see section 2.3). 

2.3. Evolution equations 

The perturbed system of a scalar field with 
a quartic repulsive interation but with temper- 
ature zero has b een studied by many authors 
(IColpi et al.ll 19861 lUrena fc BernafeOlOh . Follow- 
ing the same procedure we study the evolution of 
the SF in a FRW universe. We write the metric 
tensor as g = g° + <5g, where g° is the unperturbed 
FRW background metric, and <5g the perturbation. 
The perturbed line element in conformal time 77, 
is (we take c=l in this subsection) 

ds 2 = a(r7) 2 (-(l + 2^)dr/ 2 + 2B, i dr}dx i 

+a(T]) 2 [(1 - 24>)6 ij + 2E,ij ]dx l dx 3 , (5) 

with a the scale factor, ip the lapse function, 4> 
gravitational potential, B the shift, and E the 
anisotropic potential. We separate the Energy- 
Momemtum tensor and the field as T = To + <5T, 
and $(x M ) = $o(»7) + 5^{x^) respectively. As we 
are studying the linear regime <5$(x M ) << 3>o(»7), 
we can approximate V($>) ~ V($o). We work 
in the Newtonian gauge where the metric tensor 
g becomes diagonal and as a result, in the trace 
of the Einstein's equations the scalar potentials 
if> and <j) are identical, therefore, ip relates to the 
gravitational potential. 

Changing to the cosmological time t using the 
relation (d/dr))=a(d/dt), the perturbed Einstein's 
equations SGj = 8nGSTj to first order for an scalar 
field in the Newtonian gauge (E=0=B) are 



-8irG6p$ = 6H(6 + H6) - -^V 2 <; 

GT 

8^G$ J$, 2 = 2(4> + H(f>), i , 

Z-kGSp® = 2[4> + 3H<j)+(2H + H 2 )(j)l6c) 



(6a) 
(6b) 



with ' = d/dt and H = (In a)'. In this gauge g 
is diagonal thus (f> = tp. The perturbed density 
5p$ and the perturbed pressure Sp$ are defined in 
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terms of the perturbed energy momentum tensor 
as 

6T° = -6p 9 = -(* 6$-$lil; + V,* 5$), 

ST? = -i 
a 

5 P ^ = - - y,* <y«)«. (7) 



5T 



j 



Systems ([6]) and (7J) describe the evolution of the 
scalar perturbations. To study the evolution of 
the SF perturbations we use the perturbed Klein- 
Gordon equation 



6$ 



3H5<I> 



1 



r v 2 (5$ + v,<i> $ (5$ 



2^, $0 - 4$ o = 0. 



(8) 



Equation ([8]) can be rewritten as 
d 2 V 



□<S$ 



- 460^ = 0, (9) 



where the D'Alambertian operator is defined as 



□ 



dt 2 



1 



(10) 



Essentially equation (|SJ) represents a harmonic os- 
cillator with a damping 3_ff <5$ and an extra force 
— 2<pV,$ . Equation (jSJ has oscillating solutions if 
the term — -i?\7 2 )5$> is positive. This equa- 

tion contains growing solutions if this term is nega- 
tive, that is, if is negative enough. From here 
we see an important feature. These perturbations 
grow only if V has a maximum, even if this is a 
local one. Here the potential is unstable and dur- 
ing the time when the scalar field remains in the 
maximum, the scalar field fluctuations grow until 
they reach an stable point. This implies that the 
galactic halos could have formed within this period 
and with similar features. Finally for the back- 
ground field equation we have that 0saO, therefore 
its equation reads 



$0 + 3ff$ + V,» ($0) = 0. 



(11) 



We have that $0 depends only on time, we now 
suppose that the temperature is sufficiently small 
so that the interaction between the SF and the rest 
of matter has decoupled, and we assume that the 
symmetry break (SB) took place in the radiation 



dominated era in a flat universe, in this case the 
equation for the SF perturbations reads 

□<5$ + ^[/c|(T 2 -T^) + 12$^$_4$ o 

+ \[kl{T 2 -T 2 )+A<S>l\<S>^ = (12) 

In Figure []] we show the behavior of the poten- 
tial for different temperatures, we see that the SB 
takes place at T=Tc- 

In Magana et al. (2012b) they solve numerically 
the evolution of the fields $0 an d 8$ that satisfy 
©, ©, (HQ) and ||I2J). They show that as the tem- 
perature decreases and goes below Tc, $0 falls 
rapidly to a new minimum where it will remain 
oscillating. In a similar way, the SF fluctuation 
grows very quickly only when $0 approches the 
minimum. Taking when T>Tc, the behav- 
ior of the field $ just after the SB is what we had 
expected, i.e, it changed from being a local min- 
imum to an unstable local maximum of the po- 
tential rolling down to the minimum $ m i n . It is 
in this period when the SF fluctuation can grow. 
Hence, it is only before the SB when the SF re- 
mains homogeneous. 

From here we see that in the SFDM model, the 
primordial DM halos form almost at the same time 
due to the phase transition produced by the SB. In 
the non-linear regime these halos can merge with 
other halos forming larger structures just like in 
the CDM model. However, the main difference 
lies in the initial formation of the DM halos, they 
are formed very rapidly and almost at the same 
time, from here we expect that they possess sim- 
ilar features. This difference between the SFDM 
and CDM models can be tested by observing well 
formed high-redshift galaxies and also by compar- 
ing characteristic parameters of several DM domi- 
nated systems, for instance, by observing that in- 
deed, dwarf or LSB galaxies possess cores even at 
high-redshift, esp ecially since CDM simul ations of 
dwarf galaxies by iGovernato et al.1 (|2012h suggest 
that their DM density profiles were initially cuspy 
but later on turn into core profiles due to feedback 
processes. 

2.4. Newtonian limit 

In this work we are interested in the galaxy af- 
ter its formation. Thus, we constrain ourselves 
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to solve the static limit of equation (| 12[) when 
$ is near the minimum of the potential and af- 
ter the SB, where we expect it to be stable. We 
also expect the gravitational potential to be locally 
very homogeneous in the beginning of its collapse, 
therefore the gravitational potential <f> rj 0. The 
stability of these halos will be shown elsewhere. 
For clarity, from here on we stop using units in 
which c—1. Wc find an exact spherically sym- 
metric static solution for the SF when it is near 
the minimum of the potential (V r '($)|$ « 0), i.e., 

when $o=*min= fc K T c 7 r2 )/ 4 and T<T C . For 
the static case H 
reads 



0, $o ~ and equation ([12 



?3-(T 2 -T 2 )5<S> = 0. (13) 



Equation (fT3")) describes the evolution of the fluc- 
tuations of the SF after the SB. Moreover, this 
equation is linear. In the SFDM model these fluc- 
tuations describe the DM halos, therefore, a solu- 
tion of this system is equivalent to obtaining the 
temperature corrected density profile of early DM 
halos. 

We find that the ansatz 



~- r , sin(fcr) 
5$ = (5$ — -i — '- cos(ujt) 
kr 



(14) 



is an exaxt solution to equation (|13p provided 



; 2 2 

k c 



2ft 2 



(15) 



From equation (IT5|) we notice that now k depends 
on the temp erature T. If we use the Madelung rep- 
resentation ( Barcelo et al. 20051) we find an easier 
interpretation of the solution (|T4")) . We write 



<5$ 



2 v /gcos(Q — net) 



cos(S — fict), (16) 
(17) 



where the phases S(x,t), Q(x,t), and the func- 
tions n(x, t) and q(x, t) are taken as real functions. 
The reason of introducing k is that it gives us the 
necessary units so that we can interpret n(x,t) 
as the number density of DM particles, as $ has 
energy units. With this in mind, we can define 
an effective mass density of the SF fluctuation by 
p = mn. 



Combinig equations (fl4|) and (fTB]) and using our 
interpretation of n we obtain a finite temperature 
density profile 



P(r) = Pa 



sin 2 (fcr) 
{kr) 2 ' 



provided 



k 2 B T 2 = ft»23-4$2, 
mc 2 

S = wt + — — i. 



<1> Z 



(18) 

(19) 
(20) 
(21) 



Here k=k(T) and po=pa(T) are fitting parame- 
ters while A, Tc, k are free parameters to be con- 
strained by observations. For galaxies the New- 
tonian approximation gives a good description, 
therefore, from equation (|18[) we obtain the mass 
and rotation curve velocity profiles given by 



M(r) 
V{rf 



4xGporf sin(2fcr) \ 
~1^2V~^lkr~ (22a) 



4nGp Q 
2fc 2 



sin(2fcr) 
2kr 



(22b) 



respectively. Lately, the Einasto DM profile 
seems to give a better de scription of DM ha- 
los in CDM sim u lations ( Navarro et al. 2010t 
Merrit et all l2006t iGraham et al.l 120061 ). we in- 
clude the Einasto RC profile for a later comparison 
with ([22b]> . 



Vi =4ttGp_ 2 — 



„2/c 



a 



(3/a) 



l(-x') 
a 



(23) 

with 7 the incomplete gamma function given by 

a Jo 

and x' := l(^-) a . 

We now define the radius R of the SFDM dis- 
tribution by the condition p(R) = 0. This fixes 
the relation 



kjR = jn, j = 1,2,3,. 



(24) 



where j is the number of the exited state required 
to fit a galaxy RC up to the last measured point. 
From equations (| 13[) and (PMl) we find that the halo 
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allows exited states and then, the total density is 
the sum of densities in the different states. 

We notice that the profile in equation (Tl8|) 
presents "wiggles ," charact eristic of configurations 
in exited states (jSinl If 994). Also, we define the 
distance where the first peak (maximum) in the 



RC is reached as 



this determines the first 



local maximum of the RC velocity, which can be 
obtained form equation (j2"5)l 



(25) 



cos(27rjy) Aanfarjy) _ ^ _ Q 



2(nj) 2 y L 2njy 
where we used equation (|24|) and y := (?max/R)- 
3. Discussion 

We have seen that within our scenario of dark 
matter and galaxy formation the DM halos are 
naturally cored, i.e., their central density profiles 
are finite and do not diverge, on the contrary, they 
behave as p^r° for small r, it is important to note 
that the core is obtained from the model and not 
assumed, thus, this is an alternative way to solve 
the cusp/core discrepancy without ad ding baryons 
( Kuzio de Narav fc Kaufmann 201 lh ■ 

In the SFDM model the DM halos (density fluc- 
tuactions) form after the SB and grow only after 
the field rolls to the minimum of the potential, 
same which varies with the temperature T. There- 
fore, the initial size of the fluctuation depends on 
the temperature at which it formed. From equa- 
tion (TIB"]) we see that the size of the DM con- 
figuration, specified by R, i s now temperature- 
depen dent, therefore, as in Harko fc Madarassvl 
(|201ll ), we also solve the problem of a unique scale 
lenght for all halos, but now by using the SB mech- 
anism. Therefore, halos which formed at different 
temperatures may have different sizes (values of 
R). 

In Figure [5] we show the RC fits of three LSB 
galaxies using the minimum disk hypothesis (ne- 
glecting the baryonic com ponent) taken fro m a 
high-precision subsample of iMcGaugh I (|2005l ) , we 
compare equation (|22b|) (solid line) and equation 
(f2"3"| (dashed line) and notice that these galaxies 
present two features, long flat tails in the outer re- 
gion and wiggles. The wiggles (small oscillations) 
are perfectly reproduced by the SFDM model by 
using combinations of exited states, the letter j 
that appears in the panels of Figure 2 specifies 



the required combination of states. It is impor- 
tant to highlight the SFDM fit of NGC 1560 
we note it has the same gre at level of accurac; 



displayed in M OND models (|Gentile et al.ll201 



Sanders! [2009h . This combination of states in our 



RC fits suggests that there is not a universal DM 
profile, but instead, that the subsequent evolution 
determines the final profile. The flat outermost re- 
gion is a consequence of considering exited states, 
same behavior that was presen t in previous works 
(fSml 1 19941: iBernal et all l2010h which used T=0. 



However, the main difference now is that by con- 
sidering T^O we can accomodate exited states in 
halos and expect them to be stable due to the DM 
thermal and repulsive self interactions. Moreover, 
by considering for the first time both finite DM 
temperature corrections and exited states, we ob- 
tained an excellent agreement with the observed 
data (see Figure [2]) , suggesting that large quanti- 
ties of baryons are not essential to fit RCs in our 
model. This precludes the necessity of including 
large amounts of gas, stars, and feedback processes 
to flatten the inner regions of RCs. We expect that 
adding only the observed amount of baryons will 
be enough to reach perfect agreement with RC ob- 
servations and the correct luminosity. This will be 
seen in future works. 

Regarding the second problem about fitting the 
maximum of the RC and the outer regions at the 
same time, we notice from Figure[2]that we do not 
have such problem anymore, in contrast to previ- 
ous works (IHarko|l2011 ; Harko k. Madarassvl 201 it 
iRobles k. Matosl l2012t ). the difference is mainly 
due to the combination of states in halos. We 
can estimate a lower bound of j, the minimum ex- 
ited state necessary to agree with the data up to 
the last measured point, with the following rule 
of thumb: we find the value of y by identifying 
from the RC the last measured point as R and the 
first maximum (first visible peak) as r^ ax , we look 
for the closest value of j associated with this y in 
Figure [3] (the values of Figure [3] are determined 
by equation (|25[1 ). this provides us with the dom- 
inant state in the center of the galaxy and at the 
same time with the minimum state required to fit 
current observations. It is a lower bound because 
upcoming observations might go beyond R, in that 
case R will increase and y will decrease, implying 
larger values of j . 

As some final remarks we see from Figure [2] 
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that Einasto fits are in general agreement in the 
outer regions, but the wiggles cannot be repro- 
duced with DM only. In fact, if we want to re- 
produced the oscillations seen in high-resolution 
RCs with a non-oscillatory DM profile (NFW, 
Einasto's, Burkert's etc.), we must include the gas 
and stars dynamics in the simulations, it would be 
interesting to show the stability of the oscillations 
after including baryons, as this might be a chal- 
lenging task in LSBs galaxies due to their low gas 
content. 

In addition to the fact that the approxima- 
tion to the data by the Einasto profile is good, 
except for the wiggles and possibly the cen- 
tral region, we notice that the parameter a 
changes for each ga l axy. As noted in previous 



works flMeirit et all 120 06; 
Jin fc SutolbOOd: iGentile et al 



Navarr o et al 
20071) 



201Ct 



the change 

in a implies that halos do not possess a universal 
profile, i.e., if we take Einasto profile to be the 
best representation of the simulated DM halos in 
a CDM environment, then we should expect to 
see a non-universality in the halos of galaxies, this 
is just the same result we have obtained directly 
from the SFDM model but without assuming a 
priori a DM density profile. 

4. Conclusions 

In this paper we give an extention to the SFDM 
model that includes the DM temperature correc- 
tions to the first loop in perturbations. We pro- 
pose a Zi spontaneous symmetry break of a real 
scalar field as a new mechanism in which the early 
DM halos form. When the real scalar field rolls 
down to the minimum of the potential, the per- 
turbations of the field can form and grow. We 
give an exact analytic solution for an static SF 
configuration, which in the SFDM model repre- 
sents a DM halo. This solution naturally presents 
a flat central density profile, it can accomodate 
more than just the ground state as now the tem- 
perature T^O, and solves previous discrepancies 
in rotation curve fits at T=0, for instance, having 
a constant halo radius for all galaxies and, in the 
case of large galaxies, the incapability to fit at the 
same time the inner and outermost regions of RCs. 

Additionaly to solving these two disagreements 
we showed that there is no need to include high 
amounts of feedback to fit and reproduced the in- 



ner core and "wiggles" found in high-resolution 
RCs. Also, the model can be tested with high 
redshift observations, the SF model predicts ini- 
tial core profiles as opposed to the initially cuspy 
ones found in CDM simulation which are expected 
to flatten due to redistribution of DM by astro- 
physical processes. Finally, both the Einasto and 
SFDM RC fits suggest the non-universality of the 
DM halos, though the latter claim is still not final 
and further work has to be done in such direc- 
tion. We strongly believe that exploring further 
the SFDM looks promising to unravel the mystery 
of dark matter. 

This work was supported in part by DGAPA- 
UNAM grant IN115311 and by CONACyT Mexico 
grant 49865-E. 
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Fig. 1. — We plot the qualitative behavior of V($>) 
vs <J> for three different values of T. The solid line 
(red in the online version) is when T>Tc, here the 
system oscillates around the minimum and posses 
a Z2 symmetry, the dashed line (blue in the online 
version) is when the SB takes place, this happens 
when T=Tc, and the dotted line (colored green in 
the online version) is for T<Tc, in this period $ 
can roll down to a new minimum and there is no 
longer Z2 symmetry. 



This 2-column preprint was prepared with the AAS IATfrjX 
macros v5.2. 
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Fig. 2. — Rotation curve fits to three LSB galaxies. 
top panel: NGC 1003, middle panel: NGC 6946, 
bottom panel: NGC 1560. Solid lines (red in the 
online verison) are the fits using the SFDM model, 
dashed line (green in the online verison) represents 
Einasto's fits, and triangles are the observational 
data. In NGC 1560 we see that the dip at rss5kpc 
is reproduced more accurately in the SFDM pro- 
file. Einasto fits show different values of a in each 
galaxy, suggesting a non-universality in the DM 
halos, the same is concluded in the SFDM model. 
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Fig. 3. — We plot the relation between y and j 
obtained by solving equation (|2"5)l . Notice that for 
halos with large exited states (large j), the first 
maximum is attained at a smaller y. 
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